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ABSTRACT 

Single-stranded DNA (ssDNA) plays a major role in 
several biological processes. It is therefore of fun- 
damental interest to understand how the elastic 
response and the formation of secondary structures 
are modulated by the interplay between base pairing 
and electrostatic interactions. Here we measure 
force-extension curves (FECs) of ssDNA molecules 
in optical tweezers set up over two orders of mag- 
nitude of monovalent and divalent salt conditions, 
and obtain its elastic parameters by fitting the 
FECs to semiflexible models of polymers. For both 
monovalent and divalent salts, we find that the elec- 
trostatic contribution to the persistence length is 
proportional to the Debye screening length, 
varying as the inverse of the square root of cation 
concentration. The intrinsic persistence length is 
equal to 0.7 nm for both types of salts, and the 
effectivity of divalent cations in screening electro- 
static interactions appears to be 100-fold as 
compared with monovalent salt, in line with what 
has been recently reported for single-stranded 
RNA. Finally, we propose an analysis of the FECs 
using a model that accounts for the effective 
thickness of the filament at low salt condition and 
a simple phenomenological description that quanti- 
fies the formation of non-specific secondary 
structure at low forces. 

INTRODUCTION 

Understanding the elastic behavior of single-stranded 
DNA (ssDNA) is of great relevance because of its major 
role in many biological processes, such as replication, 
recombination, repair, transcription and transposition of 



DNA (1). In particular, characterizing the elastic response 
and the formation of secondary structure is an important 
task in comprehending the physicochemical properties 
of ssDNA. 

ssDNA forms secondary structure by Watson-Crick 
hydrogen bonding and other non-canonical base pairings 
(2). The formation of secondary structure and its elastic 
properties are also highly determined by base-stacking 
interactions (3-5). Experiments on homopolymeric 
nucleic acids have been carried out to elucidate these 
effects (6,7). The experimental techniques generally used 
for the analysis of secondary structures rely on the differ- 
ential chemical or enzymatic reactivity of single- and 
double-stranded regions of nucleic acids, but the struc- 
tural information obtained from these experiments is 
averaged over an ensemble of possible structures 
adopted by the molecule (8,9). 

Single molecule experiments (10) make possible to study 
molecular reactions by observing or measuring the mech- 
anical response of one molecule at a time. By applying 
mechanical force to the ends of an ssDNA molecule 
tethered between two surfaces, it is possible to extract 
accurate information about its elastic properties and the 
formation of secondary structure. Despite the large 
number of articles in the single molecule field, the 
existing hterature addressing the mechanical properties 
of ssDNA is scarce, especially when compared with that 
focused on double-stranded DNA (dsDNA) (11). 
Probably this is due to the fact that ssDNA is difficult 
to manipulate because it interacts with itself by forming 
hydrogen bonds (4) and binds specifically to surfaces due 
to the hydrophobicity of the exposed bases (12-15). 

Optical tweezers' experiments with ssDNA were first 
performed by force-melting a tethered dsDNA molecule 
in the presence of distilled water or formaldehyde (16) [see 
(17) for a recent overview of this method]. These experi- 
ments showed that the force-extension curve (FEC) of 
ssDNA strongly depends on the ionic strength, with 
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large deviations from an ideal chain model at high salt 
concentration. It is generally accepted that this is due to 
the formation and stabilization of secondary structure in 
high salt conditions (18). In fact, the experimental FEC at 
high divalent salt concentrations can be well described 
using theoretical models of generic semiflexible polymers 
that are capable of forming hairpins (19). To investigate 
electrostatic interactions of ssDNA, an approach that 
avoids the formation of secondary structure was imple- 
mented (18,20-22). This was achieved by using a 
chemical compound that inhibits the formation of 
hydrogen bonds (glyoxal). By screening base pairing inter- 
actions, glyoxal makes possible to study ssDNA as if it 
were a 'pure' polyelectrolyte. 

Here we use the blocking oligo method introduced in 
(23) to mechanically pull ssDNA tethers. In this method, 
ssDNA is obtained by mechanically unzipping a long 
DNA hairpin (6770 bp). Once the hairpin is fully 
unzipped and the ssDNA stretched, an ohgonucleotide 
designed to bind to the loop region is flushed in the 
fluidics chamber. The binding of the oUgo to the 
tethered ssDNA prevents the re-zipping of the hairpin as 
the tension is released. In fact, the entropic effect to close 
30 bases within a loop and the large persistence length of 
the dsDNA segment in the loop region creates a large 
kinetic barrier that must be overcome, thereby inhibiting 
the refolding of the native structure. This methodology 
allows us to study the stretching response of ssDNA 
over a wide range of forces and salt conditions and 
should be useful to study the interaction of ssDNA with 
proteins, peptides and other chemical compounds (24). 

Here we extend previous mechanical studies of the 
elastic properties of ssDNA under varying monovalent 
salt conditions (25) and report experiments at different 
divalent salt conditions. The elastic properties of ssDNA 
in presence of divalent ions have been briefly reported in 
previous works (16,18,26). To the best of our knowledge, 
this is the first work in which a systematic study based 
on force spectroscopy methods of the elastic properties 
of ssDNA at different divalent salt conditions is presented. 
Studies have been previously reported with both monova- 
lent and divalent salts in the presence of glyoxal 
(18,20-22). 

We analyze the FECs by using the most commonly 
adopted models describing the elastic properties of 
semiflexible charged polymers, such as the Worm-Uke 
Chain (WLC) (27,28) model and the Extensible Freely- 
Jointed Chain (Ex-FJC) model (16,29). These models are 
widely used to describe the elasticity of ssDNA molecules, 
but they fail to describe the FEC of ssDNA in some con- 
ditions (e.g. at low forces for low ionic strength, and at 
low and intermediate forces for high ionic strength). To 
circumvent this problem, an approach based on a polyno- 
mial expansion fit of the FEC has been previously 
proposed (30). Here we also test a description of the 
FECs based on the Thick Chain (TC) model for low salt 
concentrations and low forces (31,32). This model takes 
into account the effects of excluded volume and the elec- 
trostatic self-repulsion between the phosphate groups 
along the ssDNA backbone. Finally, we address the 
behavior of the ssDNA at high monovalent and divalent 



salt conditions with a simple phenomenological model 
accounting for the formation of non-specific secondary 
structure. 

Our work contributes to a better understanding of the 
elastic parameters describing the FECs of ssDNA at 
different monovalent and divalent salt conditions using 
analytically tractable formulas. This knowledge should 
be most useful to extract the base pairing free energies 
of dsDNA in the presence of divalent ions from mechan- 
ical unzipping experiments. This information has already 
been shown to be essential to characterize the thermo- 
dynamics of hybridization of DNA in monovalent salt 
from unzipping experiments (25). 

MATERIALS AND METHODS 

Molecular construct, experimental set up and salt 
conditions 

All the experiments were performed using a custom-built 
dual-beam force measuring optical trap as described in 
previous works (25,33). Experiments with monovalent 
ions were performed at room temperature (25° C) in a 
buffer containing TE, pH 7.5 (Tris-HCl 10 mM, 
ethylenediaminetetraacetic acid 1 mM), 0.01% NaN3 and 
a varying concentration of NaCl (10, 25, 50, 100, 250, 500 
and 1000 mM). Experiments with divalent ions were per- 
formed at 25°C in a buffer containing 10 mM Tris-HCl, 
pH 7.5, 0.01% NaN3 and a varying concentration of 
MgCl2 (0.5, 1, 2, 4, 10 mM). An additional amount of 
9mM of monovalent cations contributed by the buffer 
must be added to all NaCl concentrations and considered 
as potentially competing with divalent cations in the 
MgCh case. However, as discussed later, the latter effect 
turns out to be neghgible in the range of divalent salt 
concentrations studied in this article. The molecular con- 
struct used in the experiments consists of a long DNA 
hairpin (extracted from 1-DNA) of 6838 bp in the stem 
ending in a tetraloop. The hairpin is inserted between 
two flanking 29-bp dsDNA handles at both sides that 
provide the necessary space to carry out the manipulation 
of the hairpin [25]. By using such short molecular handles, 
the signal-to-noise ratio is increased (34). 1-DNA (New 
England Biolabs) is digested with BamHI (New England 
Biolabs) and phosphorylated at its 5'-ends. The fragment 
of the dsDNA between positions 41733 and 48502 (the 
cosR end) is used as the stem of the DNA hairpin [see 
(24) for details]. Several ohgonucleotides (sequence 
reported in Supplementary Table SI) are used as 
building blocks for the synthesis of the loop and the 
handles of the construct. The protocol for assembUng 
the loop and handles with the stem is described in the 
Supplementary Section SI). 

The two flanking handles are functionalized with biotin 
on one side and digoxigenin on the other side. 
Streptavidin-coated polystyrene microspheres (SA beads, 
3.0-3.9 i^m; G. Kisker GbR, Products for Biotechnology) 
and antidigoxigenin (Roche Applied Science) antibody- 
coated bead (AD bead, synthesized from protein 
G-coated microspheres of 3.2 ^m; Spherotech, 
Libertyville, IL, USA) are used for specific attachments 
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Figure 1. (A) Scheme illustrating the connections between a hairpin and the polystyrene beads. SA bead stands for streptavidin-coated bead, whereas 
AD bead stands for antidigoxigenin-coated bead. Beads and molecules are not in scale. (B) Steps of the protocol used in our experiments to obtain 
FECs of ssDNA. 



to the DNA molecular construct described earlier. 
Molecules are tethered to beads via specific attachments 
(biotin-streptavidin and antidigoxigenin-digoxigenin) and 
stretched using optical tweezers. SA beads are incubated 
for 20-30 min with the synthesized hairpins in TE buffer 
solution at room temperature. An AD bead is positioned 
and immobilized at the tip of the micropipette by air 
suction. Then an SA-coated bead is captured in the 
optical trap and approached to the AD bead in the 
pipette until a molecular connection is established. 
Typically, after a few attempts, a stable attachment (re- 
sisting up to a force of 40 pN) is obtained. Figure lA 
reports a schematic representation of the connections 
between the beads and the molecular construct. Next, a 
solution containing a 30-bases ohgonucleotide (blocking 
loop) complementary to the loop region of the hairpin 
(Supplementary Section SI) is flowed into the chamber. 
The oligonucleotide solution is previously incubated at 
45° C for 15 min to remove any competing secondary 
structure within the oligonucleotide. We used an oHgo 
concentration ranging from 50 to 250 nM depending on 
the salt concentration. For each ionic condition, we 
increased the concentration of the oligo from 50 nM up 
to a value for which it was possible to obtain the ssDNA 
releasing pattern in the force signal (Figure 2, red curve). 
A higher oligo concentration was often needed at low 
salt concentrations. The highest oligo concentration used 
was 250 nM for the experiments in TE 10 mM NaCl. 

Unzipping the DNA hairpin 

The procedure to measure force-distance curves (FDCs) 
of ssDNA consists of the following steps. First, a molecule 
is tethered between the bead captured in the optical trap 
and the bead immobilized in the tip of the pipette. Then 
the optically trapped SA bead is moved away from the AD 
bead at a constant puUing speed of ~30 nm/s. Meanwhile, 
the force and relative trap-pipette distance are collected 
with an acquisition frequency of 200 Hz. The tethered 
hairpin is unzipped (Figure IB, step 1) showing a charac- 
teristic sawtooth pattern in the FDC with gentle slopes 
and force rips (Figure 2, green curve). The sawtooth 
pattern is closely related to the sequence of the hairpin 
(25). The slopes are related to the elastic response of the 
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Figure 2. Cycle of pulling, showing the stretching (black) and releasing 
(red/light gray) parts of the cycle in a hairpin in TE with 100 mM of 
NaCl. The refolding of the hairpin is observed at low forces (~1 pN) as 
a sudden jump-up in the force (highlighted with an asterisk). The re- 
pulling curve in the absence of blocking oligo is drawn in green/dark 
gray [black and green/dark gray curves nearly superimpose due to the 
quasi-reversibility of the unfolding-folding process, see [(25)]. 



ssDNA generated during the unzipping process. When a 
group of base pairs is released, the tension is relaxed 
and the force drops. For the sequence used in our experi- 
ments, the mean unzipping force lies between 10 and 
20 pN varying with salt concentration (Supplementary 
Information, Figure 2). When the hairpin is completely 
unzipped, ssDNA is obtained and the force increases 
monotonically with the distance between the beads 
(stretching of ssDNA, Figure 2, red curve). 

Transient flushing of the oligonucleotide and releasing of 
the molecule 

Once a force of ~40 pN is reached (Figure IB, step 2), a 
solution of the same buffer containing the blocking loop 
ohgonucleotide is flushed gently through the central 
channel of the fluidic chamber. The two different buffers 
can be exchanged by using a valve that regulates the buffer 
entering the chamber. The beads are then steadily 
approached at a constant speed. When the oligo binds 
to the loop region, the refolding of the hairpin is inhibited 
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due to the high kinetic barrier that must be overcome to 
refold the hairpin structure. We note that the 30 bases 
oligonucleotide can only specifically bind to the loop 
region. Therefore, it should not affect the secondary struc- 
ture formation along the ssDNA, as it is not complemen- 
tary to any other region of the hairpin. In this case, the 
releasing of the molecule shows a behavior compatible 
with a molecule of ssDNA with a number of bases twice 
the number of base pairs of the hairpin (Figure IB, step 3). 
It is often possible to reach forces down to 0.5-1 pN 
without observing the refolding of the hairpin. Yet, the 
hairpin sometimes rezips before reaching such low 
forces, and a sudden jump up to the mean unzipping 
force of the hairpin is observed (Figure 2, red curve at a 
distance ~800nm, highhghted with an asterisk). By 
pulling the hairpin again, the sawtooth pattern is re- 
covered. Releasing FDCs have been collected for at least 
three molecules at each salt condition. 

From FDCs to FECs 

Conversion from distances to ssDNA molecular extension 
is needed to study the ssDNA elastic behavior. FDCs are 
converted to FECs by subtracting the elastic contribution 
of the optical trap. In the explored range of forces, the po- 
tential energy of the bead in the trap can be approximated 
by a harmonic well, so the extension or end-to-end 
distance x of the ssDNA molecule can be written as 

x^x,-(, (1) 
k 

where x,i is the measured distance, / is the force and k is 
the stiffness of the optical trap. The stiffness of the optical 
trap has been estimated for several beads at different 
salt conditions finding a value of 0.07 ± 0.005 pN/nm 
(Supplementary Section S2). Then curves have been 
ahgned using the sequence-dependent unzipping pattern 
of the hairpin previously characterized in (25) and using 
force data up to 40 pN to fit the elastic response 
(Supplementary Section S3). Finally, the FECs of 
ssDNA are ahgned with the elastic response of the 
ssDNA at the rightmost end of unzipping curves (after 
the sawtooth pattern). 

RESULTS AND DISCUSSION 

Elastic properties of ssDNA 

Results with monovalent salt ( NaCl ) 

Figure 3 shows the FECs measured at different monova- 
lent salt concentrations filtered at 2-Hz bandwidth. To 
illustrate the variability in our measurements, we show 
two different FECs at each salt condition (differences 
can arise for instance from variabihty in the geometry of 
the attachment of each molecule to the bead surface). The 
force monotonically increases with the extension. At low 
forces (<10 pN) and high salt (>100mM), a plateau com- 
patible with the formation of secondary structure is 
observed (18,19). At low forces (<5 pN), the higher the 
salt concentration, the shorter the molecular extension. 
This behavior can be attributed to the formation of 
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Figure 3. FEC of 13-kb ssDNA at varying NaCl concentration taken 
at constant pulling speed of ~30nm/s and filtered at 2-Hz bandwidth. 
Inset shows how the formation of secondary structure reduces the 
apparent contour length of the molecule (dashed lines). 
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Figure 4. FEC of 13-kb ssDNA at varying MgCl2 concentration 

taken at constant pulling speed of ~30nm/s and filtered at 2-Hz 
bandwidth. 



secondary structure that reduces the effective length of 
the ssDNA (Figure 3, inset). The increase in the height 
of the plateau at higher salt concentration shows that 
base pairing is stabihzed with salt [35]. Moreover, the 
Debye screening length becomes shorter with increased 
salt concentration (~ 0.3 nm at 1 M NaCl), and electro- 
static repulsion between the phosphate groups is 
too short-ranged to prevent base pair formation. The 
experimental FECs in Figures 3 and 4 do not show the 
characteristic features observed for single-stranded 
homopolymeric nucleic acids (6,7) that form hehces 
stabilized by the stacking of consecutive bases. This is 
not surprising if we consider that our construct contains 
a random and equally assorted amount of different 
nucleotides. 

Results with divalent salt (MgCl2) 

In Figure 4, we show the FECs at different divalent salt 
concentrations. For every concentration, two FECs are 
reported. For monovalent salt, the plateau observed 
at low forces for high salt concentrations is due to the 
formation of secondary structure. At equal molar 
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concentrations, the relative height of the plateau observed 
in the case of divalent salt is larger compared with the case 
of monovalent sah. At 10 mM MgCb, the height of the 
plateau is comparable with that observed for 1 M NaCl. 
In general, with MgCU, the plateau extends down to 
low extensions (<4|.mi), falhng to 0 force more steeply 
at low extensions. This effect is neat at high divalent salt 
conditions (Figure 4); in contrast, in the other extreme 
conditions, the equivalent FECs at low monovalent salt 
conditions for NaCl (e.g. 10 mM) hardly exhibit a plateau 
(Figure 3). 



Table 1. Persistence and contour lengths obtained from a fit of the 
FECs (in the range 10-40 pN) to the WLC model at varying NaCl 
concentrations 



NaCl 
(mM) 


(nm) 


Standard 
deviation (nm) 


(nm/base) 


Standard deviation 
(nm/base) 


10 


1.19 


0.06 


0.68 


0.02 


25 


1.04 


0.06 


0.69 


0.02 


50 


0.99 


0.05 


0.68 


0.02 


100 


0.87 


0.04 


0.69 


0.03 


250 


0.81 


0.04 


0.69 


0.03 


500 


0.78 


0.04 


0.69 


0.02 


1000 


0.76 


0.05 


0.70 


0.02 



A fixed amount of 9mM of monovalent ion concentration from the 
buffer must be added at each condition. 



Table 2. Persistence length and contour length obtained from a fit of 
the FECs (in the range 10^0 pN) to the WLC model at varying 
MgCl2 concentrations 



MgCb 


Lp (nm) 


Standard 


L, 


Standard error 


(mM) 


deviation (nm) 


(nm/base) 


(nm/base) 


0.5 


0.93 


0.06 


0.69 


0.03 


1 


0.86 


0.04 


0.69 


0.02 


2 


0.78 


0.04 


0.70 


0.02 


4 


0.79 


0.05 


0.70 


0.02 


10 


0.75 


0.04 


0.70 


0.02 



Analysis with the WLC 

To quantify the elastic properties of ssDNA, we fitted the 
FECs to the WLC using the Marko-Siggia interpolation 
formula (28): 



ksT 
L„ 



X 



X 1 

"z;~4 



(2) 



where Lp and L, are the persistence and contour lengths, 
respectively. Even if our molecular construct contains a 
short stretch of dsDNA (88 bp) serially linked to the 
ssDNA, the contribution of these dsDNA parts to the 
measured elastic response of the ssDNA is expected to 
be negligible (<1%) due to its much higher stiffness 
(50-fold per base pair) and much shorter length. FECs 
have been fitted only between 10 and 40 pN to avoid the 
presence of secondary structure. In fact, secondary struc- 
ture formation should be suppressed at forces > 10-15 pN, 
large enough to unzip DNA (36). In Tables 1 and 2, we 
report the values of the persistence length, Lj„ and contour 
length, L,., obtained from our fits. For every salt concen- 
tration, we measured the response of ~10 different mol- 
ecules, and the value of the persistence length is the 
average of the results obtained over the three best curves 
where the effect of non-specific binding of the oligo con- 
centration was found to be neghgible (see also 
Supplementary Section S7). 

In Figures 5 and 6, experimental FECs fitted to the 
WLC model in the range 10^0 pN are shown for each 
concentration of NaCl and MgCl2, respectively. The WLC 
model describes well the curves in the range of the fit, but 
fails to describe the FEC at forces <10 pN. In particular, 
at these forces, positive force deviations from the WLC 
model are observed at concentrations that are >100mM 
for monovalent salt and in the whole range of concentra- 
tions for divalent salt. Such deviations are due to second- 
ary structure formation. Moreover, at monovalent salt 
concentrations that are <50mM, some negative devi- 
ations at forces <10 pN are also observed. These can be 
attributed to excluded volume effects (18,20,21). In fact, at 
low concentrations, the Debye length can be larger than 




Figure 5. Comparison between an experimental FEC (black) and a fit to the WLC model in the range 10^0 pN (red/gray) at different NaCl 
concentrations. Results are shown for a representative molecule. 
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Figure 6. Comparison between an experimental FEC (black) and a fit to the WLC model in the range 10-40 pN (red/gray) at different MgCl2 
concentrations. Results are shown for a representative molecule. 



the persistence length of ssDNA so electrostatic repulsion 
between different DNA segments inhibits the bending of 
the filament. 

One may ask whether bond extensibility [not included in 
the WLC formula Equation (2)] could play a role in this 
force range (10^0 pN). To answer this, we have fit the 
data with the extensible version of the WLC (28). From 
this analysis, we could not appreciate a significant im- 
provement of the quality of the fits, and the values of 
the fitting parameters changed <3% for Lp and <1% 
for L,. suggesting that at force <40 pN, the inextensible 
WLC model is sufficient to describe the observed elastic 
behavior. 

Behavior of the persistence length at different ionic 
conditions 

The persistence length, Lp, for a semi-flexible polymer can 
be written as 

Lp = Ll+Uj, (3) 

where is the intrinsic persistence length and L^j is an 
electrostatic contribution that depends on the Debye 
screening length and consequently on the inverse of the 
square root concentration of ions in solution (37^1). This 
behavior has been recently tested (42-44) with the electro- 
static term, Uj oc l/[/o«]"^^, where the exponent v can take 
different values: 2 [as predicted by the Odijk-Skolnick- 
Fixman theory (37,38)], 1 [as obtained by variational 
methods (45,46)] or even <1 (47). 

Figure 7 shows Lp as a function of the inverse of the 
squared root of salt concentration. Our results are in 
agreement with previous measurements obtained by 
chemical denaturation of a 3-kb dsDNA (Figure 7, blue 
circles), albeit the present method provides significantly 
lower standard errors. In Supplementary Section S4, we 
also compare the curves of force versus relative extension 
in both cases (25). The good agreement found at high 
forces between both measurements shows that the force 
is a homogeneous function of the relative extension, a 
unique feature of ideal elastic models. The disagreement 
observed at 1 M NaCl (Supplementary Information 



[Ion*] or 100 X [Ion**] (mM) 
1000 250 100 50 25 




Figure 7. Persistence length values for ssDNA at different salt condi- 
tions in NaCl (black squares) and MgCU (red diamonds). In blue 
circles, we show data obtained from (25). JMagnesiuin concentrations 
along the x-axis have been multiplied by a factor 100. SoHd black and 
dashed red lines correspond to the linear fits to NaCl and MgCU data. 
The parameters obtained from the fit Lp = Lp(0)+A* [/on]"""'' are: 
Lp{0) = 0.68 ± 0.04 and A = 2.19 ± 0.10 (x^ = 0.666) in the case of 
NaCl; and L,,(0) = 0.70 ± 0.04 and A = 1.59 ± 0.57 (x" = 0.701) in the 
case of MgCb. 



Figure S3, right) shows that secondary structure forma- 
tion starts to be important at forces as high as 15 pN. 

The comparison between monovalent and divalent salts 
shown in Figure 7 shows a linear dependence for both 
types of salts that collapse into a single hnear trend 
if a factor of 100 between both ions is considered 
(FEC[Mgci,] ~ FECioo[Naci])- These results confirm the hy- 
pothesis that magnesium is ~ 100 times more efficient than 
sodium in non-specifically screening DNA, which is 
reflected in the elastic properties of ssDNA. Moreover, 
this effect hinders the presence of monovalent ions due 
to the protonation of Tris in the solutions with magne- 
sium. As has been suggested in the literature (22), such 
factor between monovalent and divalent salt concentra- 
tions cannot be explained by a Debye-Hiickel picture of 
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a diffuse cloud of screening counterions. However using 
strongly correlated liquid models of counterion/macroions 
interaction (48) a factor ranging from 50 to 200 is expected 
(22). Our results are consistent with findings reported in 
the case of the folding of a ribozyme (49) and in pulhng 
experiments of RNA hairpins (43,44). 

Figure 7 also shows the results of a fit with v = 1. For 
both salts, a value of Lj] equal to 0.7 nm is obtained. This 
value is in good agreement with the previously reported 
values for single-stranded nucleic acids, in particular for 
poly-U (0.67 nm) (14,50) and for chemically denatured 
ssDNA (0.6 nm) (20,21). If we perform a fit with v, a 
fitting parameter, we then obtain v = 1.15 ± 0.25 in the 
case of NaCl and v = 1.65 ± 0.85 in the case of MgCl2. 



Table 3. Kuhn length L/^ and stretch modulus K obtained by fitting 
FECs in the range 10-40 pN at varying NaCl concentrations 



NaCl (mM) (nm) 


Standard 


>f (pN) 


Standard 




deviation (nm) 




deviation (pN) 


10 1.91 


0.04 


760 


50 


25 1.78 


0.06 


670 


50 


50 1.68 


0.05 


700 


50 


100 1.55 


0.04 


630 


70 


250 1.43 


0.06 


707 


60 


500 1.38 


0.05 


710 


50 


1000 1.31 


0.06 


850 


100 


A fixed amount of 9 niM of monovalent ion 


concentration contributed 


by the buffer must be added at each condition. 




Table 4. Kuhn length Lk 


and stretch modulus K obtained by fitting 


FECs in the range 10^0 pN at varying MgCl2 concentrations 


MgCb (mM) Lk (nm) 


Standard 


A-(pN) 


Standard 




deviation (nm) 




deviation (pN) 


0.5 1.59 


0.06 


710 


60 


1 1.50 


0.03 


700 


50 


2 1.41 


0.03 


630 


60 


4 1.41 


0.05 


670 


50 


10 1.34 


0.04 


710 


60 



Analysis with the Ex-FJC 

Another commonly used model to describe the FEC of 
ssDNA is the Ex-FJC model (29): 



coth 



'/Lk 



(4) 



where is the Kuhn length, the dimension of the 
minimal segment on which the polymer can be mapped, 
and K is the stretching modulus. Similarly to (16), fits of 
the FECs were performed by fixing the contour length of 
the molecule to the crystallographic length / of ssDNA 
as described in (18) (/ = 0.57iini or L,. = 7857 nm for a 
13650 bases molecule). The fits were carried out in the 
range 10^0 pN (Supplementary Section S5), and the 
results are reported in Tables 3 and 4 and summarized 
in Figure 8. The Kuhn length decreases as the salt concen- 
tration increases with a logarithmic dependence on ionic 
strength. Our values for the Kuhn length agree with those 
reported in the Hterature, which fall in the range 
1.5-1.6nm at [NaCl] = 150mM (16,51,52) and in phos- 
phate-buffered sahne (53). Let us note that, for a fixed 
value of the contour length of the ssDNA, the introduc- 
tion of extensibihty in Equation (4) through the stretching 
modulus K implies to use a slightly shorter length of the 
crystallographic value (0.57 nm) with respect to the experi- 
mentally reported range (0.60-0.70 nm). K shows a non- 
linear trend with concentration, and ranges between 650 
and 850 pN. In the literature, values of 800 pN at 
[NaCl] = 150mM (16,51,52) and 300 pN in phosphate- 
buffered sahne (where monovalent salt concentration is 
~150mM) (53) have been reported. The values for 
given in Tables 3 and 4 and those of the persistence 
length Lp shown in Tables 1 and 2 approximately follow 
the relation Lk — 2Lp as expected for semiflexible polymer 
models when the value of is comparable with the 
segment length /. 

Low salt conditions: indications on effective thickness 

To analyze the effect of the excluded volume, FECs at low 
monovalent salt concentration (10-50 mM of NaCl, where 
secondary structure formation is not observed) have been 
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Figure 8. Kuhn length and stretching modulus values for ssDNA versus salt concentration (normal-log scale) for NaCl (black squares) and MgCl2 
(red diamonds). Magnesium concentrations have been multiplied by a factor 100. 
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10 mM 25 mM 50 mM 




Extension (nm) 



Figure 9. Fit of FECs to the TC model at low monovalent salt condi- 
tions. Experimental data (black) and fits to the TC model (red/gray) for 
low NaCl concentrations. 



Table 5. Geometric parameters obtained by fitting the FECs at low 
monovalent salt concentration with the TC model 



NaCl 


A 


St.Dv. 


L(nm) 


St.Dv. 


L, 


St.Dv. 


(mM) 


(nm) 


(nm) 




(nm) 


(|im) 


\xm 


10 


0.83 


0.03 


1.32 


0.07 


8.48 


0.10 


25 


0.74 


0.04 


1.21 


0.06 


8.57 


0.10 


50 


0.71 


0.02 


1.27 


0.03 


8.46 


0.10 



A fixed amount of 9 mM of monovalent salt concentration contributed 
by the buffer must be added at each condition. 



fitted to the TC model (31,32,42) in the whole range of 
forces studied. This model takes into account excluded 
volume effects for a polymer, and it can be used to infer 
the effective thickness of polyelectrolytes from stretching 
measurements. In the TC model, the polymer is described 
as a tube of uniform cross-section A, with ininimal 
segments of length /, persistence length Lp and contour 
length L,.. It provides a relation between the force 
applied / and the extension of the polymer x. This 
model has been used to describe the FEC of a 
homopolymeric nucleic acid Poly(U) (32). Figure 9 
shows fits to the FEC, while Table 5 reports the values 
of the fits. The values obtained are again coinpatible with 
those reported for homopolymers in the saine salt concen- 
tration range [Poly-U, A in the range 0.8-0.9 nm (32)]. 
Our results show that the effective thickness decreases 
with salt concentration owing to the more effective 
screening of the phosphate chain by cations. 

Secondary structure formation 

The blocking ohgonucleotide used to inhibit the fonnation 
of the hairpin does not prevent the formation of non- 
native secondary structure at low enough forces driven 
by stacking and base pairing interactions. In fact, the 
shoulder observed in the FECs at high salt concentrations 
suggests that ssDNA collapses in a condensed phase. This 
effect is more prominent in the presence of divalent 
cations indicating the successful coordinated action of 
the doubly charged cations to bring distal segments of 
ssDNA close to each other. In fact, stretching/releasing 
curves of ssDNA in the presence of MgCl2 exhibit larger 
hysteresis effects than for NaCl (data not shown), 
associated to the formation and rupture of secondary 
structure. 



Despite the fact that we do not have a phenomenological 
model describing this condensed phase, the large number 
of structural motifs that can be formed and their low ener- 
getic stabihty produce a soft shoulder in the FEC without 
any visible force rips. The competition between the large 
number of possible configurations that can be formed in 
such condensed phase suggests that configurational 
entropy (rather than enthalpy) is the main driving force 
toward the formation of secondary structure that 
prevents the reformation of the native hairpin. Several 
effects contribute to the soft shoulder observed in the 
FEC: the large flexibility of the ssDNA due to its low per- 
sistence length, the so-called compensation effect (54,55) 
between different structural motifs that unfold/fold inde- 
pendently due to thermal forces and the low cooperativity 
of the unfolding/folding transition expected for such 
motifs. Next, we quantify the amount of secondary struc- 
ture formed as a function of force over different salt 
conditions. 

Estimation of the fraction of unpaired bases at high salt 
condition 

To have a quantitative measure of the fraction of ssDNA 
that forms secondary structure, we used the WLC model 
to estimate the effective contour length Ljr'^' as a function 
of force. For each FEC, we extracted the effective contour 
length by fixing the value of the persistence length to 
that obtained by fitting the elastic response in the range 
10^0 pN. If the WLC model were to describe well the 
whole range of measured forces, the effective contour 
length would be equal to the one measured in the force 
range 10^0 pN (L'^^"*") and should not depend on force 
even at low forces. The fraction of unpaired ssDNA bases 
is then simply given by the ratio (p(f) = Lf^^if)/Ll^~^^, and 
results are shown in Figure 10. 

Below 50niM NaCl, L^^^ increases at forces <5 pN and 
the fraction (/>(/) becomes >1 (data not shown). This odd 
result indicates the inadequacy of the WLC model in this 
regime and is related to excluded volume effects that 
become prominent at low salt concentration, as discussed 
in the previous section. Above 100 niM NaCl, and in aU 
explored concentrations of divalent salt, Lj^" decreases 
with decreasing force. We attribute this feature to the for- 
mation of secondary structure. At low forces (<5 pN), the 
higher the salt concentration, the shorter the effective 
contour length, meaning that at high monovalent salt con- 
centration (or in presence of divalent salt), the fraction 
of paired bases increases. At the lowest measured forces 
(1-2 pN), the fraction of unpaired bases can reach up to 
10-2()% at the highest salt concentrations (namely, 1 M of 
NaCl and lOmM of MgCy. 

To extract the fraction of unpaired bases at 0 force, 
0(0), we fit data to the following phenomenological 
formula, 

l+exp[-(^)_ 

where /* is the critical force at which half of the nucleo- 
tides are paired, while 8 is a measure of the width of the 
force region in which ssDNA starts to form secondary 
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Table 6. Parameters obtained from a fit to Equation (5) of the 
curves presented in Figure 10 



NaCl (mM) 


/ (pN) 


S (pN) 


0(0) 


100 


0.440 


2.090 


0.448 


250 


3.6 


1.75 


0.313 


500 


4.850 


1.637 


0.049 


1000 


5.322 


1.520 


0.029 


MgCl2 (mM) 


/ (pN) 


S (pN) 


m 


0.5 


0.625 


1.958 


0.421 


1 


2.270 


2.070 


0.250 


2 


3.310 


1.580 


0.110 


4 


4.899 


1.610 


0.046 


10 


5.681 


1.430 


0.018 



For the NaCl case a fixed amount of 9 mM of salt concentration from 
the buffer must be added at each condition. 



Structure. An analogous exponential dependence of the 
fraction of secondary structures has been suggested on 
the basis of mesoscopic models that take into account 
the formation of multiple and degenerate structures (56). 
Figure 10 shows fits of Equation (5) to the experimental 
data, and fitting parameters are summarized in Table 6. 

For both monovalent and divalent ions, /* increases, 
whereas the fraction of unpaired bases at 0 force, (p{0), 
decreases with salt concentration, as expected. Figure 11 
shows 0(0) as a function of salt concentration. At both 
IM of NaCl and lOmM of MgCla, 0(0) has almost 
reached the minimal value of 0, meaning that at 0 force 
and high salt, almost aU bases contribute to secondary 
structure formation. The dependence of 0(0) with salt con- 
centration can be described by a power law oc (c — c*y^ 
where c* stands for a critical concentration and y > 0 is an 
exponent. From the fit (Figure 1 1), one gets y 1 for both 
monovalent and divalent salts, but with a different value 
of c* (65 mM and OmM for NaCl and MgCl2, respectively, 
see Supplementary Section S6). At high c, we find/* > S, 
meaning that we can approximate 0(0) in Equation (5) by 
exp(— /'*/3) showing that /* increases logarithmically with 
c, as indicated in the results of Table 6. 
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Figure 11. Fraction of the unpaired bases at 0 force (p(f= 0) at differ- 
ent salt concentrations. Logarithmic scale is used for both axes. Black 
squares stand for NaCl and red diamonds stand for MgCl2. 
Magnesium concentrations have been multiplied by a factor 100. 
Lines are fits to the power law dependence, 0(0) a (f — c*)"'' (see text 
for details). 



Finally, we also checked the effect of the blocking oligo 
at high salt condition (Supplementary Section S7). The 
ohgo was designed to selectively attach to the loop; 
however, it also binds to other regions along the 
ssDNA. This is confirmed by the observed decrease of 
height of the force plateau when increasing the oligo 
concentration. 



CONCLUSIONS 

In this article, we investigated the elastic properties and 
the formation of secondary structure in ssDNA using the 
approach of the blocking ohgo. We systematically 
explored a range of different concentration of monovalent 
and divalent salts observing the formation of a plateau in 
the EEC of the ssDNA, which we interpret as evidence 
of secondary structure formation. At a given force, the 
amount of secondary structure (i.e. the fraction of 
paired bases) increases with salt concentration. Fitting 
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the FECs of ssDNA to the inextensible WLC, we have 
found that the persistence length depends hnearly with 
the inverse of the square root concentration of monova- 
lent and divalent salts (31,32,40), the electrostatic contri- 
bution being apparently proportional to the Debye length. 
Interestingly enough, divalent cations have a notably 
greater effect in screening electrostatic charge than mono- 
valent salt, as has been observed in previous thermal- and 
tension-induced base pair melting. The intrinsic persist- 
ence length is ~0.7nm for both types of salt, and the ef- 
fectivity of divalent cations in screening electrostatic 
interactions appears to be 100-fold as compared with 
monovalent salt, in hne with what has been recently 
reported for single-stranded RNA (43). The values of 
the Kuhn length obtained from the Ex-FJC model 
display a similar behavior. The contour length of the 
ssDNA shows a weak dependence with salt concentration, 
suggesting that the variation of the ionic force of the 
medium does not affect the conformation of the sugars 
in the force range explored; in fact, a fraction of the 
deoxyriboses could interconvert from C3-endo (inter- 
phosphate distance 0.6 nm) to C2-endo conformation 
(interphosphate distance 0.7 nm) resulting in a change in 
the contour length (2,3). The WLC is sufficient to describe 
phenomenologically the elastic behavior of ssDNA at 
forces <40 pN, but only if the contour length of the 
ssDNA is a fitting parameter to the elastic model. In 
contrast, if we fix the contour length of the ssDNA 
using the crystallographic interphosphate distance, then 
we must use an extensible model (such as the Ex-FJC) 
to accurately fit the FECs. We also reported an analysis 
of the data based on the TC model to extract the effective 
thickness of the ssDNA at low monovalent salt condition. 
Finally, we characterized the formation of secondary 
structures from the analysis of the FECs taken at high 
salt condition and proposed a phenomenological 
formula quantifying the amount of non-native secondary 
structure. 

All these behaviors rely on a trade-off between filament 
flexibility and electrostatics. Electrostatic repulsion 
between phosphate groups along the backbone stiffens 
the filament. On increasing salt concentration, the 
backbone self-repulsion is screened and hence its persist- 
ence length decreases. Moreover, the higher the salt, the 
larger the screening of the backbone, facilitating base pair 
formation. Overall, these effects contribute to the forma- 
tion of non-native secondary structure. 

As a further step, it would be interesting to apply our 
technique to study the EEC in the presence of other kind 
of cations, both monovalent and divalent, such as potas- 
sium and calcium, or even multivalent cations, such as 
hexamine cobalt compounds and polyamines. Finally, a 
deeper understanding of the kinetics and thermodynamics 
of non-specific secondary structure formation could be 
obtained by applying the present methodology to investi- 
gate ssDNA molecules of varying length and base pair 
composition. 

SUPPLEMENTARY DATA 

Supplementary Data is available at NAR Online. 
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